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Abstract: A lattice formulation of the four dimensional Wess-Zumino model that uses 
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by performing an iterative procedure in the coupling constant. The closure of the algebra, 
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1. Introduction 

Non-perturbative studies of supersymmetric theories turn out to have remarkably rich 
properties which are of great physical interest. For this reason, much effort has been 
dedicated to formulating a lattice version of supersymmetric theories. See for example [Q]- 
H and pQ, for recent reviews. While much is known analytically, the hope is that the 
lattice would provide further information and confirm the existing analytical calculations. 

The major obstacle in formulating a supersymmetric theory on the lattice arises from 
the fact that the supersymmetry algebra is actually an extension of the Poincare algebra, 
which is explicitly broken by the lattice. Indeed, in an interacting theory, translation in- 
variance is broken since the Leibniz rule is not valid for lattice derivatives CJ. Ordinary 
Poincare algebra is also broken by the lattice but the hypercubic crystal symmetry forbids 
relevant operators which could spoil the Poincare symmetry in the continuum limit. In the 
case of the super Poincare algebra, the lattice crystal group is not enough to guarantee the 
absence of supersymmetry violating operators. Without exact lattice supersymmetry one 
might hope to construct non-supersymmetric lattice theories with a supersymmetric contin- 
uum limit. This is the case of the Wilson fermion approach for the N = 1 supersymmetric 
Yang-Mills theory where the only operator which violates the N = 1 supersymmetry is a 
fermion mass term. By tunning the fermion mass to the supersymmetric limit one recovers 
supersymmetry in the continuum limit (9|, 12]. Alternatively, using domain wall fermions 



1 13] or overlap fermions 14], this fine tunning is not required. Recently, a new lattice con- 
struction for models with extended supersymmetry has been proposed |15| ]. In this case, 
the lattice preserves some supersymmetries which are enough to reduce or eliminate the 
need for fine tunning (see also fltfl ). 

In the past the lattice Wess-Zumino model has been perturbatively studied using Wil- 
son fermions ||, ||] and adding to the action a Wilson term also for the scalar fields. In the 
continuum limit this results in a cancellation of divergences between fermion and scalar 
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fields. However, scalar and fermion renormalization wave functions in general do not coin- 
cide, due to finite contributions, thus in order to restore supersymmetry in the continuum 
limit a fine tuning of the various coupling of the lattice action is needed || . For the two 
dimensional case this problem is not present [Q], at least in perturbation theory, where the 
continuum supersymmetric Ward identities are recovered in the limit of vanishing lattice 
spacing without a fine tuning . 

More recently, a lattice Wess-Zumino model has been defined in Refs. §, § using a 
general Ginsparg- Wilson operator. In this case the supersymmetric continuum limit is 
recovered without a fine-tuning also in four dimensions Q. Moreover, this formulation al- 
lows to consider Yukawa interactions which are invariant under lattice chiral transformation 



1 18], thus it appears to be suitable for chiral theories and, in particular, for supersymmetric 
gauge theories. 

In this paper we consider the four dimensional lattice Wess-Zumino model introduced 
in Refs. || [7j and show that it is actually possible to formulate the theory in such a way 
that the full action is invariant under a lattice superymmetry transformation at a fixed 
lattice spacing. The action and the transformation are written in terms of the Ginsparg- 
Wilson operator and reduce to their continuum expression in the naive continuum limit 
a — > 0. The lattice supersymmetry transformation is non-linear in the scalar fields and 
depends on the parameters m and g entering in the superpotential. We also show that the 
lattice supersymmetry transformation close the algebra, which is a necessary ingredient 
to guarantee the request of supersymmetry. We believe that the existence of this exact 
symmetry is responsible for the restoration of supersymmetry in the continuum limit, which 
has been explicitly verified in perturbation theory in the case of the scalar and fermion two- 
point functions JjJ. 

The paper is organized as follows. In Sec. || we introduce the Ginsparg- Wilson fermion 
operator and formulate the lattice Wess-Zumino action. In Sec. [3| we show how to build up 
a lattice supersymmetry transformation that is an exact symmetry of this model. In Sec. |] 
the closure of the algebra, crucial step to be satisfied in order to impose supersymmetry, 
is shown. Discussions and outlook are summarized in Sec. [5]. In Appendices A, B and C, 
some details of the calculations are presented. 

2. The Wess-Zumino model 

The Ginsparg- Wilson relation pl| 

l5 D + Z? 75 = aD l5 D (2.1) 
implies a continuum symmetry of the fermion action which may be regarded as a lattice 



form of the chiral symmetry 18 ]. As a matter of fact, the fermion lagrangian with a Yukawa 
interaction 

C = ^Dif> + gi){P + (j)P + + P_<^ ~P_)if> , (2.2) 



1 Non-perturbative effects may produce supersymmetry breaking at finite volume ]17 
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where 

P±=2(l±75), P± = -(l±75) (2.3) 

are the lattice chiral projection operators and 75 = 75(1 — aD), is invariant under the 
lattice chiral transformation 

S^p = ie^/sip , 5tp = ^756 , 5(ft = —2ie(j) . (2-4) 

By writing ip in terms of two Majorana fermions 

ip = X + iV, ( 2 -5) 

it can be seen that the interaction term in Eq. (^^) couples the two Majorana fermions and 
therefore there is a conflict between lattice chiral symmetry and the Majorana condition 
H p. This is due to the fact that the projection operators P± depend on D. Moreover, it 
has been observed that by making the following field redefinition 

il/ = (l-^D)ifi, fi = ij, (2.6) 

the Yukawa interaction becomes 

g${P + ct>P + + P^P-W (2.7) 

and the two Majorana components of ij)' decouple. Taking advantage of this property, one 
can define the four dimensional Wess-Zumino on the lattice with Majorana fermions ||. 

We start with a lagrangian defined in terms of the Ginsparg- Wilson fermions on the 
d = 4 euclidean lattice. Our analysis is valid for all operators which satisfy Eq. ([2,1|), 
however, in the following we will use the particularly simple solution given by p(| 

where 



and 



V^0(x) = -((f>(x + afi) - (p(x)) 
a 

V*0(x) = -{<f,{x) - <t>{x - afi)) (2.10) 
p a 



are the forward and backward lattice derivatives, respectively. Substituting Eq. (|2.9| ) in 
Eq. (^^) we find convenient to isolate in D the part containing the gamma matrices and 



write 



where 



D = Di+D 2 (2.11) 



Di = a{ 1 - vxfx ) ' = vimr = 7 A • (2 - 12) 
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In terms of D\ and D2 the Ginsparg- Wilson relation (|2.l| ) becomes 

2 



D( - Di = -D x . (2.13) 

The action of the 4-dimensional Wess-Zumino model on the lattice has been introduced 
in Refs. g, |?J and can be re- written using the above notation as 

Swz = Y,{ W 1 ~ ^y lD ^ - -tfDrf + Ft(l - ^D x y l F + \mxx 
[22 a 2 2 

+m(Fct> + + gx(P + (PP + + P-tfP-)x + 9{F4> 2 + (F^)) , (2.14) 



where <f> and F are scalar fields and x is a Major ana fermion which satisfies the Major ana 
condition 

X = X T C (2.15) 
and C is the charge conjugation matrix which satisfies 

C T = -C, CCt = l. (2.16) 

Moreover, our conventions are 

c^c- 1 = -( 7m ) t 

= ( 75 ) T • (2.17) 



It is easy to see that in the continuum limit, a — > 0, Eq. ( 231 ) reduces to the continuum 
Wess-Zumino action 

S = J ^x{jd + m)x + <p ] d 2 ^ + F^F + m{F^ + {F^) 

+g X (P + 4>P+ + P-4> j P-)x + g(F<f> 2 + (i^ 2 ) f )} • (2.18) 

3. The supersymmetric transformation 

If one defines the real components by 

<!>^-L(A + iB), F^^=(F-iG) (3.1) 



the Wess-Zumino action ( 2.14[ ) can be written as 

Swz = Sq + Si n t , (3.2) 

with 

So = Y,\lx(l ~ ^Dir 1 D2X ~ ~{AD\A + BD X B) 
^-^ [22 a 
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( 1 1 

Sint = J2[ 2 m * X + m ( FA + GB) + ^ 9 ^ A + ' 75jB): 



+ — g[F(A 2 -B 2 ) + 2G(AB)] 



(3.3) 



(3.4) 



The free part of the action, So, is invariant under the lattice supersymmetry transfor- 
mation 

SA = ex = X£ 
SB = -ie^x = -iX"/5£ 
8 X = -D 2 (A - ij 5 B)e — (F — i lb G)e 
5F = eD 2X 
5G = ieD 2 -y 5 x ■ 
In fact, the variation of Sq under the this transformation is 
SS = 



(3.5) 



i-^i)- 1 ^ 

2i 



D 2 (A-i l5 B)e-(F -ij 5 G)e - ~xeB x A 



+ -XlseDxB + {eD 2X ){l - ^D^F + i(eD 2l5X )(l ~ 

a 2 2 > 



By using (|2.17| ) and integrating by part 2 , this variation becomes 

2 



1 s 1 9 2 

1 - -D X y'Dl + -A 

2 a 



B 



J2{-xe (l-^D 1 )- 1 D 2 2 + -D 1 A + ixise 
z — ' L 2 a 

X L J 

-X(l - ^D 1 )~ 1 D 2 {F - i lb G)e + XD 2 e(l - fa^F + ixD 2l5 e(l - fa)'^} 
= 0, 



where we have used the Ginsparg- Wilson relation ( p,13[) , which implies 

(i - -M-^Dl 

2 a 



V 



(3.6) 



As discussed in |7j , the variation of Si nt under fl3.5| ) does not vanish because of the failure 
of the Leibniz rule at finite lattice spacing Q. 

In order to discuss the symmetry properties of the lattice Wess-Zumino model one 
possibility is to modify the action by adding irrelevant terms which make invariant the full 
action. Alternatively, one can modify the supersymmetry transformation of Eq. (|3.5| ) in 



such a way that the action (3.2) has an exact symmetry for a different from zero . Since 



For instance, for any scalar function T one has TeD^X = X-DaJ-e. 
3 A similar attempt has been proposed by Golterman and Petcher, for the 2-dimensional Wess-Zumino 
model. 
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the transformation ( |3.5[) leaves invariant the free part of the action, this modification must 
vanish for g = 0. Therefore, we introduce the following transformation 



5A = ex = xe 

SB = -i£75X = -«X75£ 

5 X = -D 2 (A - i l5 B)e - (F - i lh G)e + gRe 

5F = eD 2X 

5G = ie~D 2 >y 5 x 



(3.7) 



where R is a function to be determined by requiring that the variation of the action 
vanishes. We make the assumption that R depends on the scalar and auxiliary fields and 
their derivatives and not on x- 

The variation of the Wess-Zumino action under the transformation ( |3,7l ) is 

SS W z = Yl {^(1 " ^D 1 )- 1 D 2 Re - mx [d 2 (A - i l5 B)e + (F — i l5 G)e - gRe 

X 

+m(AeD 2 x + Fxe + iBeD 2 j 5 x - iGx^e) + -^=x{ex + lb{£lsx))x 



- V2g X (A + i^B) \D 2 (A - i To B)e + (F - i^G)e - gRe 

[{A 2 - B 2 )eD 2X + IF Axe + 2iFBx l5 e 
+2iABeD 2lbX + 2GB X e - 2iGA{x^e)^ } . 



By using the Fierz identity, terms with four fermions cancel as in the continuum. Moreover, 
g independent terms cancel out after an integration by part, and one is left with 



SS WZ = {dX [(1 - ^D^DiR + 



mR 



9 

V2 



2 x (A + i l5 B)D 2 (A-i l5 B)e 



X D 2 {A - i l5 B) 2 e\ + V2g 2 X (A + ij 5 B)Re} 



(3.8) 



The function R is determined by imposing the vanishing of 5Swz- By expanding R in 
powers of g 

R = RW + g R& + • • • (3.9) 



and imposing the symmetry condition order by order in perturbation theory, we find 



with 



= ((\-^D 1 )- 1 D 2 + m)- 1 &L 



AL = -=(2{A + i 75 B)D 2 (A - i To B) - D 2 (A - i l5 B) 2 ) 
\J2 

l —\2(AD 2 A - BD 2 B) - D 2 {A 2 - B 2 ) 



V2 



+ 2^75 



(AD 2 B + BD 2 A) - D 2 (AB)] } 



(3.10) 



(3.11) 
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To order g 2 one has 



= -V2((l - ^Di)- 1 ^ + m)- 1 ^ + z 75 B)((l - ^D 1 y 1 D 2 + m)- 1 AL , (3.12) 



and for n > 2 



= _V2((1 - l^i)" 1 ^ + m)-\A + i l5 B)R^ . 



(3.13) 



By inserting these results in Eq. (|3.9|), the function R to be used in Eq. (|3.7|) is therefore, 
the formal solution of 



[(1 - ^D 1 y 1 D 2 + m + V2g(A + ij 5 B)]R = AL . 



(3.14) 



Notice that, from the perturbative expressions ( jlTlCD and (glaft one realizes that -R — > 
for a — > 0, since AL vanishes in this limit. Indeed, AL is different from zero because of 
the breaking of the Leibniz rule for a finite lattice spacing. 



4. The algebra 

We now study the algebra associated to the lattice supersymmetry transformation (£ 
introduced in the previous section. In particular, carrying out the commutator of two 
supersymmetries we must find a transformation which is still a symmetry of the Wess- 
Zumino action, i.e. the transformations of the fields form a closed algebra, order by order 
in g. In this section, we explicitly check this fact up to order g 1 , even though the calculation 
can be generalized to any order. 

Two supersymmetry transformations on the scalar field A give 

5i5 2 A = Si{s 2 x) 

= -e 2 [D 2 {A - i^B)ei + (F - i 75 G)ei - gRsi] 



and their commutator yields 

[5 2 , 5i]A = -2eiD 2 s 2 A + g{eiRe 2 - s 2 Rei) . 
The order g 1 of the second term on the r.h.s. of (|4.1D reads 



g(e 1 R^s 2 -e 2 R (1 h l ) 
m(l-§LY) 



(4.1) 



V2ge 2 - 



D 2 {A 2 - B 2 ) - 2(AD 2 A - BD 2 B) 



(4.2) 



s m 2 (l-fI>i) + |Z)i L" 
where we used ( |3.1C ). Then, the commutator of two supersymmetries on the scalar field A 



is 



[5 2 ,Si]A = -2ei7 M £ 2 {L> 2M ^ 

m(l-fDi) 



+ 



V2m 2 (l-%D 1 ) + lD 1 



D 2 ^(A 2 - B 2 ) - 2(AD 2lM A - BD^B)] } . (4.3) 
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Similarly, the commutators of two supersymmetries on the other fields, up to terms of order 
g , are (see Appendix A for some details) 



[5 2 ,<5i]5 = -2£ l7/t £ 2 {zVB 



+ V2g 



m(l-|Z>i; 



™ 2 (l-^l) + ^l 



[ ( J 2 ,5 1 ]F = -2£i7 M £ 2 { J D 2 ^ 



^2 



V2m2(l-f J D 1 ) + fD 1 
[5 2 ,5i]G = -2e l7At e 2 {z) 2M G 



V2g 



n 2 



m2(l -§£>!) + 2D! 



LV(AB) - (AD 2/1 P + SL» 2m A)] } , (4.4) 



,D 2M (A 2 - P 2 ) - 2(AD 2lM A - BD 2fl B)] } , (4.5) 



D 2lM (AB) - (AD 2 ^B + BD 2lM A)j } (4.6) 



and 



[S2,Si]x = -^ilfi£2^D 2fl x 

g m(l - f Pi) - P 2 



P 2 (^4 - ij 5 B)^X +(A + i^B)D 2l ^ X 



D 2 [(A-ij 5 B)^ X ] 



Therefore, the general expression of these commutators is 

[fc, 6 2 }$ = t^P*(*) , $ = (A 5, P, G, X ) , 



(4.7) 



(4.8) 



where = —2e 2 j fl e 2 and P* (<£) are polynomials in $ defined as 

P*($) = D 2 ^ + 0(g) (4.9) 
where the order g 1 contributions can be read in ( |4.3D -( {4~7"|) . We have verified that the 



closure works, i.e. the action ( 2.14 ) is invariant under the transformation 

S-^ + c^P*^) (4.10) 

up to terms of order g 1 . This calculation is sketched in Appendix B. Notice that, in the 
continuum limit D 2fl — ► and the transformation ( f4.10| ) reduces to 

+ (4.11) 

since terms in fl4.3| )-( fi77| ) of order g 1 vanish due to the restoration of the Leibniz rule as 
a -> 0. 

Higher orders in g of the transformation ( 4.10| ) can be determined by using the expres- 
sion for PW given m ( |Ig) . The proof of the invar iance of the Wess-Zumino action under 
the transformation (|3.7l ) at any order in g can be similarly performed. Indeed, the closure 
of the algebra at any order in g should hold since the supersymmetry transformation (^ 
is an exact symmetry of the lattice Wess-Zumino action. 
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5. Conclusions 

In this paper, we have presented a lattice formulation of the four dimensional Wess-Zumino 
model with an exact supersymmetry using Ginsparg- Wilson fermions. We have shown 
that it is actually possible to formulate the theory in such a way that the full action is 
invariant under a lattice supersymmetry transformation at a fixed lattice spacing. This 
supersymmetry transformation introduces a function R which is non-linear in the scalar 
fields and depends on the parameters m and g entering in the interaction part of the action. 
The action and the transformation, which have been written in terms of the Ginsparg- 
Wilson fermions, reduce to their continuum expression in the limit a — > 0. We have also 
shown that the lattice supersymmetry transformations close the algebra, as it is required 
by the supersymmetry. While the present work is confined to the proof to the order g , 
concerning the closure of the algebra, there are no obstructions to extending this procedure 
to higher orders in g. 

The study of the Ward identities associated to the exact lattice supersymmetry we have 
introduced can be done by generalizing the analysis performed by Golterman and Petcher 
in ||] for the two dimensional Wess-Zumino model. In the Appendix C, we have calculated 
a simple Ward identity up to order 0(g) and verified that it is satisfied. We believe that 
the lattice supersymmetry we have introduced automatically leads to a restoration of the 
continuum supersymmetry without additional fine tuning. Explicit results on the two point 
functions jjj lend support to this idea and we are currently investigating on this issue in 
more detail. 

Obviously one the most important question is whether these ideas may be extended 
to supersymmetric gauge theories where we expect that the Ginsparg- Wilson relation will 
play an important role in the construction of a lattice supersymmetry. However, there is 
an important difference between gauge theories and the Wess-Zumino model. The free and 
the interaction terms of a lattice gauge action are both contained in the plaquette term 
and therefore it is not obvious how to perform the iterative construction of the lattice 
supersymmetry transformation in this case. 
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Appendix A 

In this Appendix we show some details concerning the calculation of the commutator 
[$2, $i]x- Applying two supersymmetry transformations on x ° ne has 

hfax = -D 2 [siX ~ 75(£75X)]^2 - [(eiAsX) + Ts^i^Tsx)]^ + g(SiR)e 2 , (5.1) 

where (SiR) is the supersymmetry variation of the function R. A similar expression is 
obtained for 525\x with e 2 £\. Terms to the order g° can be treated as in the continuum 
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case by using a Fierz rearrangement as well as the relations in Eq. ( 2.17 ), For instance, 



taking the first term in Eq. (|5.1| ) and the corresponding one with e 2 «-> £\ one has 

-S\D 2lJL x{l^2) a + £2D 2fl x( r 7^l)a = ^(£l7^2)(£>27^X)a - ^l7A^ e 2)(£>27A^X)a • (5-2) 

Using a similar rearrangement for the remaining terms order g°, the commutator of two 
supersymmetry transformations on X is 

[fa, fa]Xa = -^ilti£2D 2fl Xa + 0(g) . (5.3) 
The contribution to this commutator to the order g 1 is 

(5 1 R w e 2 - 5 2 R w ei) = -^-e 2 (fix) + -^-£2(ei75X) - (^l <-> £2) 



where the sum is over the 16 independent 4x4 matrices. By using Eq. (2.17) only the 
terms with = {7^,7^} survive, moreover, using the explicit form for one finds 
Eq. (O). 



Appendix B 

The prove of the invariance of the action under the transformation ( [4.9D to the order g° is 
immediate. In this Appendix, we explicitly calculate the variation of the fermionic part of 
the action ( |3.2[ ) to the order g 1 (the remaining part, containing only scalar fields, can be 
similarly treated). 



l ^X(l - | Al) 1 AsX + \™XX + ~^9X(A + »75#)x j 
X (1 - o-^l) D 2 + m ] 



V2^\ A v 2 i; ; m 2 (1 _| Dl) + ^ 1 

x [ D 2 (A - HsB)^ X + (A + n 5 B)D 2 ^x ~ D 2 [(A - i l5 B)^x] 



2x(A + i lb B)D 2llX ~ X(D 2 »A + i l5 D 2fl B) X , (5.5) 



where (|4.3| ), (t4.4|) and (|4.7D have been used. Due to the Ginsparg- Wilson relation ( 2.1$) , 
we have that 

and fl5.5| ) becomes 

-~J= XI {xAjO 4 - 275^)7mX + xG 4 + i^B)D 2 ^y^x ~ XI ^(A - ij5X)D 2 x 

- 2 X (A + ij 5 B)D 2lxX ~ x(D 2fl A + i-ft-D^B)*} (5.7) 
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where in the last term of the first line an integration by part and the relation ( 2.17 ) has 
been used. Finally, terms containing the derivative of the scalar fields cancel out, while to 
prove the cancellation of the remaining ones a Fierz rearrangement is needed. 



Appendix C 

By making a variation of the generating functional with respect to the lattice supersym- 
metry transformation ( |3.7[) , the Ward identity reads 



< J $ <S$ >= (5.8) 

where J$ are the sources for the fields = (A, B, F, G,x)- A simple Ward identity is 

< D 2 (A - iysB) > + <F> -i 75 < G > -g < R >= . (5.9) 

The first term of this Ward identity is zero because of 5— momentum conservation and 
D 2 (k = 0) = 0. In the following, we show that this Ward identity is satisfied to order 0(g), 
i.e. 

< F > (1) -z 7 5 < G > (1) -g < R > (0) = (5.10) 

where < O >^ denotes the expectation value of the function O to the n order in pertur- 
bation theory. From the action (|3.2|), the free propagators are 



< AA ><°) = < BB >(°)= -M-\l - -Di)- 1 

< FF >(°) = < GG > (0) = -M~ 1 D 1 

a 

< AF >(°) = < BG > (0) = mM 



-i 



< XX > (0) = -M-H(l ~ \DiY 1 D 2 - m) . (5.11) 
The term < F x >W in momentum space reads (a factor g/y/2 is omitted) 



< F(k) >W= f { 

Jpq I 



<x(p)x(q) > {0) <F(k)A(-p-q) >(°) 

pq L 

+ < F(k)F{-p - q) >(°) (< A(p)A(g) > (0) - < B{p)B{q) 
+ 2 < F(Jb)A(-p - g) > (0) < F(p)A(g) > {0) 

+ 2 < F(Jb)A(-p - q) > (0) < B{p)G{q) > (0) | . (5.12) 

By substituting the propagators we have 

< F(k) >«= ^ / ^^(pjrrCCl - ^-D 1 (p))- 1 D 2 (p) - m) 
m J p 2 

+ 4<5 4 (fc) / TW" 1 ^) = 0. (5.13) 
Jp 

The one point function of G is zero at this order. 
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Finally, the last term of the Ward identity ( |5,9[) is 



< R >W= ((1 - -D l y 1 D 2 + m)- 1 < AL > 



(0) 



(5.14) 



In the following we consider only the contribution from the field A since B can be treated 
similarly. In momentum space we have 



Indeed, the first integrand is an odd function of q, while the second term is zero since 
D 2 (A; = 0) = 0. 
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